
************************************************************
*  Regressions in Table 3
************************************************************

clear
# delimit ;
version 11;
capture log close;
set more 1;

log using table3.log, replace;
use data_1yr_stacked;

gen yob=1959-age;
gen time=yob-1915;

*  Birth rate;
*  Women age 18-44;
reg brate sr10a lnpopm urbshf i.yob i.regno, cluster(regno);
reg brate sr10a lnpopm urbshf i.yob i.regno [aweight=popf], cluster(regno);
reg brate sr10a lnpopm urbshf i.yob i.regno i.regno#c.time, cluster(regno);
reg brate sr10a lnpopm urbshf i.yob i.regno i.regno#c.time [aweight=popf], cluster(regno);

*  Urban age 18-44;
reg brateu sr10u lnpopmu i.yob i.regno, cluster(regno);
reg brateu sr10u lnpopmu i.yob i.regno [aweight=popfu], cluster(regno);
reg brateu sr10u lnpopmu i.yob i.regno i.regno#c.time, cluster(regno);
reg brateu sr10u lnpopmu i.yob i.regno i.regno#c.time [aweight=popfu], cluster(regno);

*  Rural age 18-44;
reg brater sr10r lnpopmr i.yob i.regno, cluster(regno);
reg brater sr10r lnpopmr i.yob i.regno [aweight=popfr], cluster(regno);
reg brater sr10r lnpopmr i.yob i.regno i.regno#c.time, cluster(regno);
reg brater sr10r lnpopmr i.yob i.regno i.regno#c.time [aweight=popfr], cluster(regno);

*  In logs;
reg lbrate sr10a lnpopm urbshf i.yob i.regno i.regno#c.time, cluster(regno);
reg lbrate sr10a lnpopm urbshf i.yob i.regno i.regno#c.time [aweight=popf], cluster(regno);

*  Urban age 18-44;
reg lbrateu sr10u lnpopmu i.yob i.regno i.regno#c.time, cluster(regno);
reg lbrateu sr10u lnpopmu i.yob i.regno i.regno#c.time [aweight=popfu], cluster(regno);

*  Rural age 18-44;
reg lbrater sr10r lnpopmr i.yob i.regno i.regno#c.time, cluster(regno);
reg lbrater sr10r lnpopmr i.yob i.regno i.regno#c.time [aweight=popfr], cluster(regno);

log close;


